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Abstract 

We explore the constraining power of OPE associativity in 4D Conformal Field Theory 
with a continuous global symmetry group. We give a general analysis of crossing symmetry 
constraints in the 4-point function (</></><^ '0'), where (f) is a primary scalar operator in a given 
representation R. These constraints take the form of 'vectorial sum rules' for conformal 
blocks of operators whose representations appear in R ® R and R® R. The coefficients in 
these sum rules are related to the Fierz transformation matrices for the R <8> R <8> R <8> R 
invariant tensors. We show that the number of equations is always equal to the number of 
symmetry channels to be constrained. We also analyze in detail two cases — the fundamental 
of SO(N) and the fundamental of SU(N). We derive the vectorial sum rules explicitly, and 
use them to study the dimension of the lowest singlet scalar in the <f> x $ OPE. We prove 
the existence of an upper bound on the dimension of this scalar. The bound depends on 
the conformal dimension of (j) and approaches 2 in the limit dim(0) — > 1. For several small 
groups, we compute the behavior of the bound at dim(</>) > 1. We discuss implications of 
our bound for the Conformal Technicolor scenario of electroweak symmetry breaking. 
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1 Introduction 



Conformal Field Theory was originally conceived in four and three dimensions, with applications 
to particle physics and critical phenomena in mind. However, it is in 2D that the most spectacular 
results and exact solutions have been obtained. In higher dimensions, there seems to be a general 
feeling that the constraining power of conformal symmetry by itself is insufficient to tell nontrivial 
things about dynamics. Hence the interest in various additional assumptions, like supersymmetry, 
or integrability for the planar J\f = 4 super Yang-Mills, or the AdS/CFT duality. This is not fully 
satisfactory, since there are likely many 4D CFTs which do not fulfill any of these assumptions. 
For example, "conformal windows" of non-supersymmetric gauge theories. 

And yet, in the early days of 4D CFT, it was hoped that the OPE associativity is such a 
strong constraint on the CFT data (the spectrum of operator dimensions and the 3-point function 
coefficients) that it could allow for a complete solution of the theory. Recently [1] , [2] , [3] , [I] we 
have been taking a fresh look at this idea, originally proposed by Polyakov [5]. Our approach is 
not to try to solve the OPE associativity, but rather to try to deduce from it general bounds that 
any CFT must obey. We discovered that such general bounds do exist! The bounds found so far 
fall into two general classes: 

• PQ>[2] An upper bound on the gap in the operator spectrum: any unitary J^D CFT containing 
a scalar operator of dimension must also contain another scalar O appearing in the 
OPE 0x0 whose dimension do is bounded by a universal function of d^: 



The function f{d^) has been computed numerically by means of a well-defined algorithm. 

• Upper bounds on the 3-point function coefficients X^o where is the same scalar as above 
and O is any primary operator in the OPE 0x0. They have the form: 



where g is a universal function of the dimensions and of the O's spin lo- For lo = this 
function was evaluated in [3]. Very recently, the case lo = 2, do = 4, corresponding to the 
stress tensor OPE coefficient, was considered in @,pE]. In this case the upper bound @ 
gives a lower bound on the central charge of the theory. 



do < fidf). 



(1) 




(2) 
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In this paper, we will discuss a generalization of bounds of the first class to the case when CFT 
has a continuous global symmetry G (Abelian or non-Abelian) , and an operator transforms in 
a nontrivial representation R of G. We will consider the OPE <fi x <p* if R is complex, or <fi x </> if 
R is real. We will be discussing how to show that this OPE necessarily contains a singlet scalar 
operator S whose dimension ds is bounded by a universal function which depends only on 0's 
dimension and transformation properties: 



Thus the novelty with respect to PQ,[2] is that we will be bounding the gap in a given global 
symmetry sector (singlet in this case). 

It is useful to recall that the original motivation of [1J was to find a bound of precisely this type 
for the case G = 5*0(4) and <fi in the fundamental. This in turn was needed in order to constrain 
the Conformal Technicolor scenario of electroweak symmetry breaking [7]. This connection was 
discussed extensively in pQ, and we will come back to it in the discussion Section. 

The paper is organized as follows. We first derive the "vectorial sum rules", which generalize 
the main equation encoding the OPE associativity — the sum rule of [1] — to the globally symmetric 
case. Not surprisingly, these new sum rules take a different form depending on G and R under 
consideration. To begin with, we treat three concrete examples: in the fundamental of SO(N); 
4> charged under a U(l); <fi in the fundamental of SU(N). Here we provide explicit derivations, 
illustrating the necessary technical ingredients. We then consider the general case in some detail, 
and in particular show that the number of constraints from crossing symmetry is always equal to 
the number of unknown functions. 

Then we discuss what our vectorial sum rules imply for the bounds on the singlet dimension. 
In principle, functions fs can be computed by a straightforward generalization of the algorithm 
of pp. However, numerical difficulties involved are much greater in the present case, because one 
is working in a linear space of a much bigger dimension. As a result, we cannot yet push our 
analysis to the point where it produces numerically significant bounds. Thus we follow a more 
modest strategy. First, we explain how one shows that a bound does exist for sufficiently close 
to 1. The main idea is to consider the case d^ — 1 and then to argue by continuity (which is 
possible since the vectorial sum rules are continuous in d^). This analysis also shows that fs 
approaches 2 as d^ — > 1. In other words, just like in our previous work, the free field theory limit 
is approached continuously. 

The next question is then to determine how fast fs approaches this limit. Here we limit 




(3) 
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ourselves to quoting several f$ values at d^ > 1, leaving the determination of a detailed shape of 
these bounds to future work. 

We conclude by discussing consequences for phenomenology and promising research directions. 

Note. As this work was being prepared for publication, a very interesting paper [6] appeared, 
which gave several important generalizations of our method and results. In particular, [6] gener- 
alized our method to the M = 1 SCFT case, and used it to derive bounds on non-BPS quantities. 
They also derived lower bounds on the central charge and on the two-point functions of global 
symmetry currents, with or without supersymmetry. Finally, they presented a set of equations 
incorporating OPE associativity constraints in the case of a U(l) global symmetry. They have 
even derived a bound on the U(l) singlet dimension (in the supersymmetric case). We will be 
commenting on this partial overlap in more detail below. 



2 Sum rules in CFTs with a global symmetry 
2.1 Conventions 

We begin with some preliminary comments and notational conventions. We will work in the 
Euclidean space. Just as in our previous work, availability of explicit expressions for 4D confor- 
mal blocks given by Dolan and Osborn |8J will play a crucial role. Consider a 4-point function 
(0(xi)x^(^2)x( a; 3)^( a; 4)) where <f> and x are two primary operators, not necessarily Hermitean, 
assumed to have equal dimensions d^ = d x = d. The OPE x x will contain a sequence of spin 
I, dimension A primary fields Oa,i,' 

^xx^^AaA,!. (4) 

A,/ 

Here Aa,; are the OPE coefficients, in general complex. We then normalize the conformal blocks 
via: 

4>(xi), .0^x4) \ v-^ 1 , , 

r O2) x{x3) I 12 34 

U = X 12 X 34: / (x 13 X 2 4) = ZZ, V = £14X23/ ^13^24) = (1 — z)(l — z) , (6) 

zz 

g&,l{ U i V ) = —[kA+l(z)kA-l-2{z) -{Z<r> Z)] , (7) 

z z 

k /3 (x) = x^ 2 2 F 1 {/3/2,/3/2,/3;x) . (8) 
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The points Xi are assumed to be near the vertices of a square, as the picture suggests. The 
ordering is important. Eq. ^ says that the exchanges of Oa,i and of its conformal descendants 
in the (12)(34) channel (=s-channel) can be summed up in a 'conformal block' g& t i(u,v). The 
coefficients pa,z are given by 

Pa,i = \Xa,i\ 2 > 0- (9) 

Compared to [8], and also to our previous work, we have dropped the (—1/2)' prefactor in the 
expression for g>A,z- This normalization is more convenient for the following reason. In the new 
convention all conformal blocks are positive when operators are inserted at the vertices of a square 
in the shown order (this corresponds to z — z — 1/2). This is just as it should be, because this 
configuration is reflection-positive in the Osterwalder-Shrader sense with respect to the vertical 
median line (notice that the fields in the two sides of the correlator are complex-conjugate of each 
other )Q Thus any s-channel contribution to the correlator, even spin or odd, has to be positive. 
There is no disagreement with Doland and Osborn |Sj, because in their notation the extra minus 
sign would be offset by a change in the sign of the OPE coefficient in the RHS of the correlator. 

The (14) (23) channel (=t-channel) conformal block decomposition can be analyzed similarly. 
In this case we will need OPEs x 0^ and \ x x and only fields appearing in both of these OPEs 
will give a nonzero contribution, proportional to the product of the two OPE coefficients. 

In p] we have analyzed the particular case when is Hermitean and \ — 111 this case the 
s- and t-channels correspond to the same OPE (0 x 0). In addition, only even spins contribute 
because of permutation symmetry x\ -H- X2- Let us introduce the notation for the sum of all 
s-channel contributions: 

G+ = Pa,i9a,i(u,v) , (10) 

l even; A 

(+ means that we are summing over even spins only) and a tilde notation for a contribution of 
the same set of operators in the t-channel: 

G + = G+^ v = ^2pa,i9a,i(v, u) . (11) 

Here we used the fact that going from the s- to the t-channel, which means simply rotating the 
picture by 90°, interchanges u and v. In this notation the crossing symmetry constraint of p] is 
written compactly as: 

G + ={^ d G + . (12) 

1 Actually, conformal blocks are positive on the whole interval < z = z < 1. Configurations corresponding to 
such z, z can be mapped onto a rectangle, which is reflection-positive. 
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The appearance of the (u/v) d factor in this relation is due to a nontrivial transformation of the 
prefactor 1 / (x^x^) in ([5j under crossing. 

Finally, an important technical remark. Unlike in pQ, to extract full information from the 
4-point function ([5]), we will have to consider not only the s- and t-channel OPEs, but also 
the u-channel ones (13) (24). Conformal blocks for such 'diagonal' OPEs are related to the 
nearest-neighbor conformal blocks discussed above by analytic continuation, which introduces 
spin-dependent signs into the crossing-symmetry constraints. A useful way to keep track of these 
signs is not to consider the u-channel OPE directly, but to instead apply the s- and t-channel 
decompositions to the 4-point function with the permuted insertion points: 



Here, we transposed the fields in the right side of the correlator. Now in the t-channel we have the 
same OPE as we would have in the u-channel in ([5]). And in the s-channel we have the same OPE 
as in ([5]), except for the transposition. This transposition is taken into account by reversing the 
sign of the odd-spin contributions in the s-channel (and permuting the flavor indices accordingly, 
see below). 



As a first example we will now consider the SO(N) global symmetry case, with a scalar primary 
operator <p a transforming in the fundamental representation. We normalize the 2-point function 
of (p a as (<p a (x)(pb(0)} = 5 ao (x 2 ) d , d = d^. Consider the 4-point function 



The operator insertion points are assumed numbered in the same order as in ([5]). 

Operators appearing in the a x fa OPE can transform under the global symmetry as singlets 
S, symmetric traceless tensors T( afe ), or antisymmetric tensors A[ ab y. 




(13) 



2.2 Fundamental of SO(N) 




(14) 



(j) a X (j) h = 5 ab t 



(15) 



+ SabS^ (even spins) 
+ T^l (even spins) 

+ A \ab\ (° dd S P inS ) • 
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The index (a) shows that an arbitrary number of operators of each type may in general be present, 
of various dimensions A and spins I. However, permutation symmetry of the a 0& state implies 
that the spins of the S"s and T's will be even, while they will be odd for the A's. 

We note in passing that the stress tensor will be an S of A = 4, I = 2, while the conserved 
SO(N) current will be an A of A = 3,/ = 1. The OPE coefficients of these operators are related 
to the stress tensor and the current central charges by the Ward identities [11] . which should allow 
to derive various bounds on these quantities by the method of [3]. The simplest cases of these 
bounds have already been explored in [I], [6]. In this paper we will not be making any assumptions 
about the central charges of the theory and will treat the stress tensor and the current on equal 
footing with all the other fields. (However, in future studies central charge information may be 
useful; see the discussion Section.) 

On the other hand, it will be important for us that the unit operator 1 is always present in 
the (f) a x 4>b OPE, with a unit coefficient. 

As mentioned in the Introduction, we are interested to learn something about the dimension of 
the lowest-dimension singlet scalar (an S of I = 0). This will require disentangling its contribution 
to the 4-point function from the possibly present low-dimension scalars of type T. 

We will now see what the crossing symmetry says about the relative weights of various con- 
tributions in the 0x0 OPE. Applying the conformal block decomposition in the s-channel we 
get: 



a d 
b c 



I I'< i+G *>+(„ + X4l I)-g t+ (_-x)' g -< i6 > 



Here Gs,t,a are defined as in (10), and sum up conformal blocks of all fields of a given symmetry. 
Remember that Gs,t contain only even spins, while Ga only the odd ones. The unit operator 
contributes together with the singlets, and its conformal block is = 1. To keep track of the index 
structure, we are using the graphical notation for tensors. Every line means that the corresponding 
indices are contracted with the 5 tensor: 

J [ = 5 ab 5 cd , etc. (17) 

The index structure of the symmetric traceless and the antisymmetric tensor contributions in 
(16) is fixed by the symmetry (and by the tracelessness, in the case of Gt)- The signs are fixed 
from the requirement that for a = d ^ b = c all contributions have to be positive by reflection 



positivity, see Section 2.1. Apart from the sign and the index structure, we do not keep track of 
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the overall, positive, normalization of each term. In other words, we know that each G contains 
conformal blocks summed with positive coefficients, but we do not keep track of the normalization 
of these coefficients. This is sufficient for deriving constraints on the operator spectrum, which is 
the focus of this paper. On the other hand, normalization conventions will be important for any 
future studies of the OPE coefficients. 

Next we apply the t-channel conformal block decomposition to the same 4-point function, and 
we get an alternative representation: 



a d 
Ibc J 



uy 
v J 



[l + Gs)+([ \ + 



X 



2 
N 



I I 



X 



G, 



Note that to get this equation requires only changing the index structure appropriately, permuting 
u 4-¥ v (here we are using the tilde notation introduced in Section 2.1), and multiplying by {u/v) d 
to take into account how the 1 / (a^f 2^34) transforms. 



Now we equate the s- and t-channel representations and pick up coefficients before each of the 

. We get 2 independent equations: 



3 inequivalent tensor structures: , * | 



u 



u- d {G T -G A } 



1 + Gs — j^Gt 



,-d 



Gt — Ga I 

Gt + Ga 



and a third one which can be obtained from the second by u -H- v : 

2 



1 + Gc 



N 



G T \ = u - d {Gt + G A } 



(18a) 
(18b) 

(19) 



Notice that for the SO(N) case using the u-channel OPE would not yield any new equation. 



It will be convenient to rewrite the system (18a), (18b) in the following equivalent form: 

. 0, 

2 ' 



F s + 1 



He 



1 + 



N 
2 

N 



F T — Fa 
F t + F A 

H T — H A 



(20a) 
(20b) 

(20c) 



where we introduced notation for (anti) symmetric linear combinations of G and G: 



Thus (20a) is obtained from ( 18a) just by grouping and dividing by if 



-u 



. Eq. (20b) is obtained 



by taking the difference of (18b) and (19) and moving the contribution of the unit operator to 



the RHS. Finally, Eq. (20c) follows by taking the sum of (18b) and (19), and again separating the 
unity contribution. 

Note that the functions F(u, v) were already used in our previous work, while the appearance 
of H(u, v) is a new feature of the global symmetry analysis. Writing the equations in terms of 
these functions is convenient because they are highly symmetric with respect to the z — z — 1/2 
point (they have only even derivatives in z + z and z — z at this point). 



The system (20a)- (20c) is then the main result of this Section. In an expanded notation, it 



can be written 'vectorial sum rule': 



ZX 



/ \ 

Fa,i 
Ha,i 



( 



+ XX 



( -F AJ \ 



XX' 



v 



—Ha l 



J 



( \ 
1 

-1 



(22) 



v-v 



Fa,i 

\-(l + %)H Ajl ] 

Here the functions F& t i(u,v) and H& t i(u,v) are related to the individual conformal blocks g^i by 
the same formulas as F and H are related to G. Their dependence on d is left implicit. In each 
sum we are summing vector-functions corresponding to the dimensions and spins present in this 
symmetry channel, with positive coefficients. The total must converge to the constant vector in 
the RHS. 

Consequences of this new sum rule for the lowest singlet dimension will be discussed below. 
Let us do however a quick counting of degrees of freedom. In total we have three G-functions: 
Gs,Gt,Ga, each of which is restricted only to the odd or even spins. The vectorial sum rule gives 
three equations for their (anti) symmetric combinations F and H. This coincidence between the 



number of equations and unknowns is not accidental; see Section |2.5| One may hope that the 
constraining power is similar to the case without global symmetry, when we had one equation for 
only one function G + . We will see in Section^ how this hope is realized. 



2.3 U(l) 

We next discuss the U(l) global symmetry, intermediate between SO(N) and SU(N). On 

the one hand, we will be able to check that the U(l) constraints agree with the already considered 
SO(N) case for N = 2. On the other hand, the derivation will be similar to the SU(N) case which 
follows. In particular, we will be working with complex fields and will need the u-channel OPE. 
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We want to derive constraints from crossing in the 4-point function of a charge 1 complex 
scalar 0. Charge normalization is unimportant. The nonvanishing correlators must have zero total 
charge, thus we are led to consider (cjxjxft . There are two basic OPEs: 



Charge sector: x $ = 1 + spins 0,1,2 . . 
Charge 2 sector: 0x0 = even spins only . 



(23) 
(24) 



Let us begin by considering the configuration 



(25) 



which is the same as in ^ for x = 0. By doing the s- and t-channel conformal block decomposi- 
tions and demanding that the answers agree we get a constraint: 



u 



{1 + G+ + G } =v~ d {l + G+ + G j . (26a) 
Here the subscript refers to the charge fields appearing in the relevant 0x0^ OPE. As indicated 



in (23), this OPE contains both even and odd spin fields, whose contributions we separate in Gq. 



According to the discussion in Section 2.1, reflection positivity of (25) implies that even and odd 



spins contribute in (26a) with the same positive sign. 



Next consider the configuration with the transposed right side of the correlator: 



Equating the s- and t-channel decompositions we get: 



u 



{1 + G+- Go} 



(26b) 



The LHS of this equation differs from the LHS of (26a) only by the reversed sign of the odd spin 



contribution (see Section 2.1). The t-channel decomposition appearing in the RHS is positive 



since the configuration is reflection-positive in this direction. 



Eqs. (26a), (26b) solve the problem of expressing crossing constraints in a U(l) symmetric 
theory. Very recently, the same equations also appeared in [6]. The authors of [6] have noticed 



that they could get a bound on the lowest dimension singlet by using just Eq. (26a) (they only 
computed the bound in the supersymmetric case, but the general case must be similar). Dropping 
the other equation simplified the problem, but the downside was that they had to Taylor-expand 
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up to a pretty high order (12) to extract the bound. Below we will show that if one uses all 
equations the second-order expansion is already sufficient to extract a bound. 

Upon identification 

Gs = Gq, Ga = G , G T = -G2 (27) 

the U(l) constraints become equivalent to the N = 2 case of the SO(N) constraints discussed 
above. The appearance of a positive factor 1/2 is consistent with the fact that we are keeping 
careful track of positivity but not of the normalization. 

2.4 Fundamental of SU(N) 

Our last example is the SU (N) case, with a primary scalar (pi transforming in the fundamental. 
We have two basic OPEs: 

(pi x </>! = Sal + S iT x Singlets(spins 0,1,2. . . ) + Adjoints(spins 0,1,2. . . ) , (28) 

's (odd spins) . (29) 



<pi x (pj = 's (even spins) + 



The representation content of the first OPE is N ®N = 1+Adj. Notice that, in general, there will 
be singlets and adjoints of any spin. The adjoint sector will contain the conserved current, but at 
present we are not using information about its coefficient. The second OPE contains symmetric 
and antisymmetric tensors, of even and odd spins respectively. 

The constraints are now derived by a combination of what we did for SO(N) and U(l). First 
consider the following 4-point function configuration: 

0i. ,4>\ 

4' '<P 3 

The s- and t-channel conformal block decompositions are evaluated using the first OPE. Equating 
them, we get a constraint: 



u~ d 



1 ^)+ I ^ -^1 )(Gt^ g- X( ,: 



= v- d {*^(i+G: s +G s ) + (\ I-^^)( G i dj +G Adj )} 



Here lines denote SU (iV)-invariant contractions of iV (dots) and iV (circles) indices by S^. The 
tensor structure of the Adj contributions is fixed by the tracelessness condition of the SU(N) 
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generators. The sign is fixed by the condition that for % = j ^ j = % the s-channel contributions 
must be positive by reflection positivity. 



Setting equal the coefficients before 



and 



we get two equations: 



u 



1 



l + G+ + G 5 --(G+ dj + G 



Adj, 



G Adj + G 



Adj f > 



and a second one which is just the u ■<-)■ v version of the first. 
Next we consider the transposed 4-point configuration: 



(30a) 



Equating the s- and t-channel decompositions, we get: 



-d 



G~s) 



X 



1 

iVi 



(^Adj G Ady 



X 



Gr 



The s-channel decomposition is obtained from the previous case by transposing the index structure 
and flipping the sign of the odd-spin contributions. The t-channel decomposition is obtained by 



using the second OPE (29). The index structure is fixed by (anti) symmetry of the exchanged 
fields, while the signs are determined by demanding positive contributions for i — t ^ j — J 
(which makes the configuration reflection-positive in the t-channel). 

Collecting coefficients before the two inequivalent tensor structures, we get two more equations, 
which this time are independent: 



u 



u 



1 + G + S -G- S 



{^Adj ^Adj} 



N G tdj + ^Adj 



v - d i G m - a 



m 



Gr 



Bj ' 



(30b) 
(30c) 



The system (30a)-(30c) solves the problem of expressing the crossing symmetry constraints. 



Like in the SO(N) case, we will find it convenient to rewrite it by separating the unit operator 
contributions and (anti) symmetrizing with respect to u and v. We end up with the following 



11 



equivalent vectorial sum rule: 



F s + F s +(l-i^A dj +(l-^)^A dj = 1 

H + s +H S + -(1 + ^)^ = "I 

F s ~ F s ~ii F Adi +jj F Adj + F tn = 1 

-^Adj -*g = o 

^Adj _iJ Adj -% = °- 

Just like for SO(N), the number of components, six, is equal to the number of the OPE channels 
classified by representation x (spin parity): S ± , Adj ± , m, []. 



2.5 General case 

In this Section we will consider the case of an arbitrary global symmetry group G, with (f> a 
transforming in an irreducible representation R. We aim at a general analysis of crossing symmetry 
constraints. In particular, we would like to understand why the number of constraints came out 
equal to the number of unknown functions in the explicit SO(N) and SU(N) examples above. 

We will assume that R is complex. The case of R real is analogous but simpler; necessary 
changes will be indicated below. 

To understand the group theory aspect of the problem, we begin by counting the number 
of scalar invariants which can be made out of two 0's and two 0^'s. These invariants can be 
constructed by decomposing the products (j) a x </4 and (f)p x </>t into irreducible representations 
and contracting those. The tensor product representation decomposes as: 

n 

R®R = ($r i (+f i ), (32) 
i=i 

where (+fj) indicates that the representations in the RHS must be either real or come in complex 
conjugate pairs. To simplify the discussion, assume for now that all r\ are real and different. In 
accord with the above decomposition, we have 

^x4 = EE^n. (33) 

i Ai 

where the objects \& a, transform in the r i: and C l aBlA . are the Clebsch-Gordan coefficients (A { 
is the index in the r,). Then we can construct exactly n invariant tensors by contracting two 
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Clebsch-Gordan coefficients: 

so that the product of two 0's and two 0*'s can be decomposed into a sum of T's: 

j 

where in the second line we indicated that we can do the same construction in a crossed fashion, 



invariant tensors satisfy a linear relation ('Fierz identity') 



by starting with the <f) a x 0^ product. The fact that both decompositions exist means that the 



Taa/3/3 ~ -^i'^ap^a " ( 3 ^) 

The matrix T is invertible and must satisfy T 2 = 1, since crossing is a Z 2 operation. 

It is also possible to construct invariants by starting from <p a x 0^, which requires the tensor 
product 

n 

R®R = Q)fj. (37) 

3=1 

Assume for now that all f's appearing in this product are also distinct (excluding as well the 
possibility for the same representation to occur both in the symmetric and antisymmetric part 
of the tensor product). Under this simplifying assumption, the number of f's is the same as the 
number of r's. Indeed, we can construct invariant tensors 

where C J aj3A , (resp. C^ A ) are the Clebsch-Gordan coefficients for fj in R x R (resp. fj in R x R). 
These must be related to T's by another Fierz identity 

where T is again an invertible matrix. Notice however that T and thus T 2 ^ 1. 

After this prelude, we come back to our problem of analyzing the crossing symmetry constraints 
of the CFT 4-point function. 
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Step 1. Let us compare the s- and t-channel conformal block decompositions: 

ft 8 




E 




(40) 



ft 



P 



Introduce functions Gi which sum up conformal blocks of operators in the representation r% (which 
will in general occur in both even and odd spins). The tensor structure of these contributions will 
be given precisely by the invariant tensors T introduced above. The crossing symmetry constraint 
then takes the form: 

E T aappGi(u, v) = J2 T* am Gi(v, u) . (41) 

i i 

Here we assume that the signs of T's have been chosen in agreement with reflection positivity. To 



simplify the notation we included the u 
separate the unit operator explicitly. 



prefactors in the definition of Gi. We also do not 



Eq. (41) will be consistent with the first Fierz identity (36) if and only if 

Gi(u,v) = Tl'Gi>(v,u) . 
Let us now define even and odd combinations: 

( ~ ±) G t = G l (u,v)±G l (v,u). 



(42) 



(43) 



These are the analogues of the F and H functions from Eq. (21). We put the index (±) on the 
left to stress that it has nothing to do with the spin parity index used in the previous Sections; 
these functions receive contributions from both even and odd spins. We have 







(44) 



where P± = (1 =F J 7 )/2 are projectors, (P±) 2 = P± by using T 2 = 1. Going to the diagonal basis 



for P±, it is clear that Eq. (44) represents a total of n constraints. 

Step 2. We next compare the s- and t-channel conformal block decompositions of the transposed 
4-point function: 



a 



E 




ft 



p 




(45) 
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The crossing symmetry constraint can be written in terms of the invariant tensors introduced 
above as: 

i 3 

Here we have shown explicitly that the odd spin parts G~ of the Gi flip signs compared to the 
above configuration (40). Note as well that each of the functions Gj will include even or odd spins 



only, depending if r, occurs in the symmetric or antisymmetric part of R x R. 



For Eq. (46) to be consistent with the second Fierz identity (39), we must have 



Gf{u, v) - Gr( u , v) = T/G^v, u) . (47) 

Since the functions in the RHS and LHS now refer to completely different OPE channels (r* in 
(j) x <fy vs fj in x 0), this equation gives exactly 2n constraints when (anti) symmetrizing in u, v. 

To summarize, we expect 3n constraints for 3n channels rf ,fj. In particular, n = 3 for the 
fundamental of SU(N). 

In case when R is a real representation, we only have one set of invariant tensors, whose Fierz 
dictionary matrix satisfies T 2 = 1. In this case each of n representations in the R x R product 
will contribute with only even or odd spins. Only the first step of the above analysis is needed in 
this case. We will get n constraints for n channels. The fundamental of SO(N) corresponds to 
n — 3. 

Generalizations . Let us now discuss how one can relax the assumptions on the content of R®R 
and R(g) R taken in the above argument. In general, R® R may contain repetitions of the same 
representation as well as conjugate pairs, while R ® R may contain the same representation in 
both symmetric (s) and antisymmetric (a) part. As it will become clear below, these two things 
must happen simultaneously. A sufficiently representative example is R = 15 of G = SU (3) [9] : 

15 ® 15 = 1 + 64 + (8i + 8 2 ) + (27j + 27 2 ) + (10 + TO) + (35 + 35) , (48) 
15 (8) 15 = 3 a + 6 S + 15' s + 24 a + 42 a + 60 s + (15, + 15 a ) + (21 a + 24 s ) . (49) 

In 15 ® 15 we have 8 and 27 appearing twice each, and also two conjugate pairs (10 + 10 and 
35 + 35), while in 15 £g> 15, 15 and 24 appear both as s and a. 

In cases like this, it is slightly more involved to count the quartic invariants. When counting 
in the R<S> R channel, every conjugate pair r + f gives two invariants which for future purposes 
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we (anti) symmetrize with respect to (aa) <-> ((3/3): 

^aaA^ppA ± ^aaA^p^A ' l 0U J 

A 

In the same channel, a fc-fold repetition of a real representation r gives rise to k 2 invariants: 

E^w^ (i,j = i...fc), (51) 

which can be (anti) symmetrized with respect to (aa) -H- producing A; (A; + l)/2 symmetries 

and A;(A; — l)/2 antisymmetries. 

When counting in the R® R channel, every representation r occurring both as s and a gives 
rise to 4 invariants 

out of which two are symmetric and two antisymmetric in (aa) -H- (/3/3). 

The total number of invariants must of course be the same counted in R <g> .R and in R® R 
channel. This is indeed true in the above example, when both 15 £g> 15 and 15 ® 15 give 14. The 
number of symmetric in (aa) -H- (/3/3) invariants also agrees (10 in both channels). This is also 
true in general. An intuitive argument is as follows. The total number of invariants equals the 
number of independent coupling constants in the scalar potential V((pi,fa 2 ,<p3,fa) where fa are 
four non-identical scalars transforming in R. This number should be the same whether you begin 
by contracting fa with <p\ or 03. Analogously, the number of symmetric invariants is the number 
of quartic couplings if we identify fa = fa, fa = fa. 



Each of the two Fierz identities (36) and (39) will now split into two, one for symmetric and 
one for antisymmetric invariants. 



Let us now proceed to the crossing symmetry analysis of the 4-point function y , ; , , 
To begin with, out of all the invariant tensors discussed above, only the symmetric ones will 
appear as the coefficients in the conformal block expansions of this correlatoij^] The (aa) 4-> ((3(3) 
symmetry is made manifest by applying a conformal transformation which maps a generic 4-point 



configuration in (40) onto a parallelogram. The 180 rotation symmetry of the parallelogram then 



acts on the indices as (aa) ((3(3), see Fig. [T] 



2 In a general Lorentz-invariant theory, the flavor structure of this correlator will involve both symmetric and 
antisymmetric tensors. 
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-T <fh& 4 ) 
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<t>U X 2) 1 ' M X 3> 

Figure 1: For any 4-point configuration, there exists a conformal transformation which maps it 
onto a parallelogram. 

To see how this symmetry arises in the conformal block decomposition, consider the OPE 

K x 4 = £ ( E ^C^a) O a (53) 

r real \ i=l / 

+ E x oc: aA o A + (-iyx* c: aA o\. (54) 

r+f pairs 

Here in the first line we include all operators belonging to the real representations. If the rep- 
resentation is repeated k times in the R® R product, there are k independent Clebsch-Gordan 
coefficients, and k independent real OPE coefficients \ % . In the second line we have operators 
from the complex-conjugate pairs, whose OPE coefficients are always complex-conjugate, up to a 
spin-dependent minus sign. 



By using this OPE in the s-channel conformal block decomposition of (40), we see that indeed 
only symmetric invariant tensors arise. Notice that the off-diagonal invariant tensors (i ^ j) in 
the case of repeated representations will appear with coefficients AqAq (x conformal block), which 
are not positive definite. We will discuss below what this means for the subsequent application of 
the derived constraints. 



We then consider the t-channel conformal block decomposition of (40), and repeat the analysis 
of Step 1. The resulting number of constraints is equal to the number n sym of symmetric invariants, 
while the number of representation x (spin parity) channels is 2n sym . 



To generalize Step 2, we have to consider the t-channel decomposition of (45). In this channel, 
the OPE parity selection rules imply immediately that only symmetric tensor structures appear, 
in agreement with the above general result. If the s and a representations are not repeated, as in 
the 15 ® 15 example, then only diagonal terms are present, and all conformal blocks enter with 
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positive coefficients. Compared to Step 1, we will have n sym new representation x (spin parity) 
channels and 2n sym new constraints. 

We are done: we have a total 3n sym constraints for 3n sym channels. Moreover, these constraints 
distinguish not only different representations appearing in the OPE, but also different copies of 
the same representation, and how they 'interfere' among each other. 

Let us now come back to the fact that if repeated representations are present in R ® R, the 
off-diagonal 'interference' channels have coefficients AqAq. To appreciate the difficulty that this 
creates, readers unfamiliar with our method of linear functionals are encouraged to read the rest 
of this Section after having read Section |3j 



Consider then our abstract way (58) of representing the vectorial sum rule. It is crucial for 
us that when all coefficients p a are allowed to vary subject to the positivity constraints p a > 0, 
linear combinations in the LHS fill a convex cone. In particular, this allows us to use the dual 



formulation of the problem in the form (75). Since the off-diagonal coefficients may be negative, 
the geometric interpretation in this case is not as obvious. Notice however that the off-diagonal 
coefficients cannot become arbitrarily negative since they are not independent of the diagonal 
ones. For a sharp formulation, consider a symmetric real matrix 

o 

where we allow for presence of more than one operator O with a given dimension, spin, and 
representation. The characterizing property of P is positive-definiteness: 

PijSiSj^O Vsj 6 R . (56) 

Now, as can be seen from this equation, the set of positive-definite matrices forms by itself a convex 
cone. It follows that the set of vectors in the LHS of the vectorial sum rule will remain a convex 



cone even if repeated representations are present. Constraints (56) replace the simple inequality 
Pa > 0. In practical applications these constraints may have to be discretized by choosing a finite 
set of vectors s,-. 



The dual formulation (75) is extended to the present case as follows. For the vectors Xjj in 
the LHS of the sum rule corresponding to diagonal (i = j) and off-diagonal (i ^ j) channels of 
the repeated representation, the simple condition A[x a ] > must be replaced by the following 
condition on the matrix A[xy]: 

PijA[xij] > VP positive-definite. (57) 
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In other words, A[Xy] must belong to the cone dual to the cone of positive-definite matrices. 
However, the latter cone is in fact self-dual, as can be easily inferred from the representation (55). 
Thus, A[xy] must be itself positive-definite. 

In the above discussion, only real representations were allowed to repeat in R <g> R. However, 
repetitions of complex pairs could be treated similarly; the only difference is that the corresponding 
P matrices will be positive-definite Hermitean rather than real. 



3 Bounds on the lowest singlet scalar dimension 
3.1 Generalities 

The previous Section would be a futile exercise in group theory if our vectorial sum rules did not 
have any useful consequences. We will now discuss how they can be used to bound the gap in the 
singlet scalar sector. Consider the SO(N) case for definiteness. Given a CFT spectrum, the sum 



S T A 

rule (22) can be viewed as an equation for the coefficients p A y . If we start imposing restrictions 



on the spectrum, such as raising the singlet scalar gap, it is conceivable that this equation will 
not have any solution consistent with the positivity requirement p&,i > 0. This is in fact precisely 
what will happen. 



3 



An equivalent, geometric, way to view this is as follows. Let us rewrite the sum rule (22) 
schematically as an equation in a linear space V of functions from two variables u, v into R 
(vector space of vector-functions): 

^p Q x Q = y. (5 



Here vectors x a represent all vector- functions appearing in the LHS of (22), while the y is the 
vector corresponding to the RHS. 



For a fixed CFT spectrum and varying p a > 0, the vectors in the LHS of (58) fill a convex 
cone, so the question is whether the vector y belongs to this cone. Imposing restrictions on the 
spectrum reduces the set of vectors x Q generating the cone, and the cone shrinks. It may be that 
the new smaller cone no longer contains y, see Fig. [2]^a,b). 

In this paper we will be dealing with two types of restrictions on the spectrum. First of all, 
we will always impose the unitarity bounds [10] 



A > 1 (1 = 0), A > I + 2 (I > 1) . (59) 
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(a) (b) (c) 



Figure 2: Geometric interpretation of the sum rule: (a) the sum rule has a solution y belongs 
to the cone; (b ) the assumed spectrum is such that the sum rule does not allow for a solution 
y does not belong to the cone; (c) for Ag = A", the y belongs to the cone boundary. 

These lower bounds on operator dimension are a completely general property of unitary 4D CFTs. 
Note that they depend only on spin and not, say, on the global symmetry representation in which 
the operator transforms. 

Second, we will impose a lower bound on the dimension of scalar singlets: 

A > A 5 (/ = singlets only). (60) 

According to the above discussion, increasing As makes the cone shrink. Our goal will be to show 
that for Ag above a certain critical value A c s r , the y is not in the cone. This critical value will 
then be a theoretical upper bound on the dimension of the first singlet scalar, valid in an arbitrary 
unitary 4D CFT. This is the bound denoted by fs in Eq. ^ of the introduction. 

Note that the list of vectors generating the cone, and thus the cone itself, vary continuously 
with As. Since for A s < A 5 r the y is inside the cone and for A s > A c s r it's outside, for A s = A s r 
it must belong to the cone boundary, see Fig. [2^c). 

Up to now we were keeping the external scalar dimension d<p fixed. However, the vectors 



entering the sum rule depend on d^ via Eq. (21). This dependence is continuous, and thus the 



cone will vary continuously with d^. It follows that the bound fsid^,), if it exists, will have a 
continuous dependence on d^. 

The reader may be worried that the above discussion was not totally rigorous. Indeed, the 
vector space V is infinite-dimensional, and there may be subtleties of convergence. However, below 
we will always be considering a finite-dimensional subspace of V, by Taylor-expanding the sum 
rule up to a fixed finite order k. On the one hand, this means that the bounds that we will derive 
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will not be optimal (they will approach optimality in the k — > oo limit). On the other hand, 
finite-dimensional analysis is simpler both practically and from the point of view of mathematical 
rigor. In particular, the statement that fs is a continuous function of d<f> is safe at finite k. 



3.2 Best possible bound for = 1 

We will begin by analyzing carefully the case d<f, = 1. The reader may wonder why this is 
necessary, since there is a theorem that a scalar saturating the unitarity bound is necessarily free. 
This theorem is easy to prove: starting from the 2-point function (0(x)0(O)) = x~ 2 we deduce 
(d 2 (fi(x)d 2 0(0)) = and thus d 2 <p = in the Hilbert space. In the SO(N) case, we will have N 
free real scalars <p a . The <p a x 4>b OPE will contain two scalar operators 

= :0 C C :, = :<Pah - (61) 

the first a singlet, the second a symmetric traceless, both dimension 2. Bounding the singlet scalar 
dimension from above seems a waste of time: the exact solution shows that it must be exactly 2. 

Yet it is not useless to give a derivation of such a result — that at = 1 there must be a singlet 
scalar of dimension < 2 — using our method based on the vectorial sum rule. The reason is very 
simple: unlike the derivation sketched above, our method is robust with respect to small variations 
in df/). If we show that it works at d^ — 1, it is guaranteed to work as well for d^ sufficiently close 
to 1. Thus we will know that a bound exists for close to 1 and, since it is a continuous function 
of d^, that it approaches 2 as d^ — > 1. We believe that such an existence proof is conceptually 
important. It is also easier than actually computing the bound at d^ > 1. The latter problem will 
be discussed in the next Section. 



Consider then the sum rule (22) for = 1. This sum rule has at least one solution-the one 
corresponding to the theory of N free scalars. This solution is very special, in that among all the 
vectors appearing in the RHS, only those corresponding to twist A — I = 2 fields will have nonzero 
coefficients. This is because in free theory no fields of other twists appear in the <fr a x <p b OPE. 
Apart from the two A = 2 scalars mentioned above, there are infinitely many twist 2 fields of the 
form 

■4a d M1 . . . d m (j) b : , (62) 
appropriately (anti) symmetrized in a, b to separate SO(N) representations. Expanding the free 
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field theory 4-point function 



S S 

((p a (xi)4> b (x 2 )(l)c(x 3 )(f) d (x^) = I \ +crossings (63) 



/-v>2 ,-y»2 

x 12 x 34 



into twist 2 conformal blocks, one can find all the p&,i coefficients. The S and T contributions 



can be disentangled since they have a different index structure, see Eq. (16). Even though we will 
not need the explicit expressions, we quote the result: 

Pl + 2,i = (^)7(2/)!, pf +v = (2/N)pJ +2;l (I even) , 

p*. v = (l\) 2 /W (/odd). (64) 

Notice that all p^,i > 0, consistent with unitarity. 

We now proceed to showing that any solution of the sum rule at = 1 must contain a singlet 
scalar of dimension < 2. In fact we will show an even stronger result — that any such solution 
must contain a singlet scalar of dimension exactly 2. 

Let's group the twist 2 terms in (p2|) separately from the rest: 



twist 2 twist ^2 

First we will show that in any solution all twist ^ 2 coefficients pp must be zero. This is shown 
by exhibiting a linear functional Ao : V — > R such that 

A [y] = o, 

A [xJ = V fields of twist 2, (66) 
Aofc/?] > V fields of twist ^ 2. 



Applying A to (65), and using the fact that pp > 0, we conclude immediately that all pp = 0. 

The functional Ao can be written explicitly by Taylor-expanding the vectors entering the sum 
rule around the point z — z — 1/2. Let's introduce the coordinates a, b 

z = i + a + b, z = ^ + a - b . (67) 

Then the functions F^j and H^,i are even with respect to both a and b, so that their derivatives 
d™ d]!; at a = b = are nonzero only if both m and n are even. It turns out that the functional A 
can be chosen as the following linear combination of second derivatives at the a = b = point: 

A [x] = A(d 2 aXl - d 2 bXl ) + B(d 2 a x 2 - d 2 b x 2 ) + C(d 2 a x 3 - d 2 b x 3 ) . (68) 
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Here £1,2,3 are the components of the vector-function: x = (xi, X2, x%) T . There is a certain freedom 
in choosing the coefficients A,B,C; the following choice is one possibility which works for all 
N >2: 

A = l, B = 2, C = 0. (69) 



The first relation (66) is trivially satisfied; the other two follow from the following curious property 



of the .Fa,; second derivatives at the a = b = point: 

%F Ail = dlF^i at twist 2 , 

dlF A}l > d 2 b F^ at twist ^ 2 . (70) 

Since all pp = 0, we are reduced to a simpler equation which involves only twist 2 fields: 

Po^ a = y • (71) 

twist 2 

We want to show that if we drop the scalar singlet from this equation, there are no solutions. This 
is shown by exhibiting a second linear functional A x with the following properties: 

Ai[y] <o, 

Ai[x a ] > on all twist 2 fields except for the singlet scalar. (72) 
Notice that any such Ai must necessarily be negative on the singlet scalar, to allow at the very 



least the existence of one explicit solution (64). 



This functional can be written again as a linear combination of derivatives at a = b = 0: 

AiM = Yl E §iW*<- (73) 

A — 4 *■ — 4 mini 

i = l m,n even 
0<m+n<k 

At present we can only find the coefficients A^ n numerically. For 2 < < 7 it is enough to 
use derivatives up to the second order, just like in A . For example, for A^ = 4 one can use the 
functional whose only nonzero coefficients are 



^0,0 - 


-8, 


A 0,2 


= 6, 


A o,o - 


-11, 


X 2 


= 8, 


X 3 - 
A o,o - 


-10, 


X 3 
A 0,2 


= 3. 



(74) 



Including derivatives up to the fourth order (k = 4) allows to find functionals in the range up 
to A^ < 128. While we have not checked higher N, we feel sufficiently confident that, adding 
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more and more derivatives, functional Ax can be found for any N. With this small proviso, the 
demonstration that at — 1 there must be a dimension 2 singlet scalar is complete. 

It is perhaps useful to give a geometric representation of the given proof, see Fig. [3j The 
existence of Ao means that there is a hyperplane (the zero set of Ao) such that all the twist 2 
vectors as well as the vector y belong to it, while all twist ^ 2 vectors lie on one side of it. The 
existence of Ai means that this hyperplane can be rotated so that the twist 2 singlet scalar and 
the rest of the twist 2 vectors lie on the opposite sides of the rotated hyperplane. 

TQ = 0) 




Figure 3: This figure gives a geometric interpretation of the proof that at — 1 the sum rule 
has no solution unless the A = 2 singlet scalar is included in the spectrum. The solid-contour 
plane represents the zero set of the functional Ao. The vector y and all the twist 2 vectors (black 
dots) lie in the A = plane. On the other hand, all twist ^ 2 vectors, which for varying A trace 
separate curves labeled by spin and representation, lie on one side of this plane > 0). The 
Ao = plane can be slightly rotated so that the y vector and the twist 2 singlet scalar lie on one 
side of the rotated plane, while the rest of the twist 2 vectors lie on the opposite side. The rotated 
plane (dashed contour) can be described by an equation A + eAi = for a small e. 



Let us now discuss the SU(N) case. The statement is the same: at d^ — 1 the OPE <fi x 



must contain a dimension 2 singlet scalar. This is derived from the SU(N) vectorial sum rule (31) 
by using the same method of A and Ax. Just like for SO(N), the functional A can be given by 
using only second derivatives. For Ai, one can use derivatives up to the second order for A" = 2,3, 
while derivatives up to the fourth order work for at least all A^ < 50. 

Finally, the reader may want to compare the above discussion with Section 5.4 of PQ, where 
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the d<f, = 1 case was analyzed analogously, and somewhat more explicitly, for CFTs without global 
symmetry. 



3.3 Some numerical results at gL > 1 



As explained in Section 3.1, the bound f$ is a continuous function of d^. The previous Section 
has established that fs = 2 at d^ = 1. The next question is to understand how fast the bound 
deviates from 2 at we increase d^ above 1. To do this we need to be able to compute the bound at 
any given d^. This is done by the method of linear functionals pQ, already used above to analyze 
the d<f, = 1 case. At any fixed d^, we will be looking for a linear functional A such that 

A[y] < 0, 

A[x Q ] > V scalar singlets with A > As and (75) 
V other fields (subject to the unitarity bounds) . 
If such a functional exists, then the spectrum without any singlet scalars below As cannot solve 



the vectorial sum rule and cannot be realized in any CFT. Geometrically, Eqs. (75) says that there 



is a hyperplane separating the y vector from the cone generated by such a spectrum. The bound 



fs is computed as the smallest As for which a functional satisfying (75) can be found 



Assuming that A is of the form (73), Eqs. (75) define a linear programming problem for the 



coefficients \ % mn - However, the number of constraints is formally infinite, which requires careful 
discretizations and truncations. We refer to PQ,[2], an d also to jB], for a detailed description of 
how these numerical difficulties are overcome. 



There is no difference of principle between the system (75) and the analogous system analyzed 
in the case without global symmetry in our previous work; it is just bigger. There are several 
families of vectors x a , labeled by representation and spin parity, and each vector now represents a 



vector-function. As a result, for the same number of derivatives (parameter k in (73)), the globally 
symmetric case involves Q times more constraints and Q times more functional coefficients, where 
Q = 3; 6 for SO(N); SU(N). This makes computations much more time-consuming as well as 
increases numerical instabilities. In the case without global symmetry, the initial analysis of [1] 
used k — 6, and subsequently we were able to push k up to 18 [2J, producing a very strong bound. 
For the globally symmetric case, we have so far not been able to go beyond k — 4 for SO(N) and 
k = 2 for SU(N). 

Here is an account of these exploratory calculations: 
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• In every case that we looked at, we found a bound for near 1 which approached 2 
continuously in the — > 1 limit. Thus we have checked the existence theorem from the 
previous Section. 

• We have seen that the bound is a monotonically increasing function of at d^ > 1. 

• Typically, we have only computed the bound in the interval 1 < d^ < d* where fs(d*) ~ 4. 
This is because knowing whether the singlet scalar is relevant or irrelevant is a particularly 
interesting question. It is also sufficient to get an idea about how strong the bound is, for 
a given k. In the following table, we give the for the cases that we considered. At the 
current level of accuracy, our f s interpolates almost linearly between 2 and 4 as d^ increases 
from 1 to d*. 



G 


U(l) = SO{2) 


SO(3) 


SO(4) 


SU(2) 


SU(3) 




1.063 (k = 2) 
1.12 (k = 4) 


1.032 (k = 2) 
1.08 (k = 4) 


1.017 (k = 2) 
1.06 (k = 4) 


1.016 

(k = 2) 


1.003 

(k = 2) 



Table V. First line: the global symmetry groups considered. The external scalar was assumed to 
transform in the fundamental. Second line: the value of external scalar 's dimension d^ for which 
the bound fs on the singlet scalar dimension (monotonically increasing from 2 for d^ = 1) was 
seen to cross 4. 



As we already mentioned in Section 2.3, a U(l) bound on the singlet scalar dimension was 
recently published in [6]. It is not possible to compare our and their results directly: on the 
one hand, they assume supersymmetry and expand to a very high order (k = 12) which should 
make their bound stronger, but on the other hand they use but one scalar component of the full 
vectorial sum rule, which makes their bound weaker. It seems that these two effects compensate 
each other, so that their bound is roughly comparable to our U(l), k = 4 bound. Notice however 
that they could not see any bound for d > 1.16, while we checked by using the full vectorial sum 
rule that a bound continues to exist even for larger d. 



4 Discussion 

In this paper we extended the constraints from OPE associativity and crossing first derived in 
Ref. PQ to CFTs with global symmetry. Focussing on the scalar 4-point function we derived 
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a set of sum rules that constrains the operator content in all the possible channels with given 
symmetry and parity (— l) e quantum numbers. More precisely, by explicit examples and by a 
general argument, we have shown that the number of sum rules equals the number of possible 
channels. This results suggest that, in principle, one could obtain (correlated) constraints on the 
operators appearing in each channel. In analogy with previous studies these constraints could 
involve the dimension and the fusion coefficients of the lowest lying operators in a given channel. 

As a first exploration we have studied the possibility of obtaining an upper bound on the 
dimension of the scalar singlet of lowest dimension appearing in the OPE of <p x <fi and <px<p\ That 
question is also relevant to asses the viability of the so called Conformal Technicolor (CT) scenario 
[7] . The goal of CT is to achieve a natural separation between the electroweak scale and the scale 
of Flavor. That is motivated by the experimental success of the Cabibbo-Kobayashi-Maskawa 
pattern of flavor violation. The way CT works is that the electroweak symmetry breaking sector 
above the weak scale flows to a strongly coupled CFT. The role of the Higgs field is then played 
by a composite operator H, while electroweak gauge interactions arise from the weak gauging of a 
SU (2) x U(l) subgroup of the global symmetry group G of the CFT. In that situation, suppression 
of flavor violation is the more robust the closer to 1 is the dimension du of H. On the other hand, 
a natural separation of mass scale requires all total singlet scalars to be only marginally relevant 
or irrelevant. Indicating by S the lowest dimension singlet C H x W, we thus need ds slightly 
below or above 4 |^J A major constraint on du is imposed by the top coupling, that runs like 
At(/i) = \t{nEw){nl ' HewY h ~ 1 and quickly becomes strong for dn — 1 = 0(1). Notice that one 
major difficulty to achieve the dream of CT is that for 4 = 1™ must have ds = 2 and the 
hierarchy problem in all its splendor. 

To be more quantitative about the needed pattern of field dimensions, assumptions on the 
physics of flavor must be made. Making the optimistic, although plausible, assumption that flavor 
violation in the light families is either suppressed by their mixing to the third family of by their 
Yukawa couplings^} the range dn ^ 1-7, ds <^ 3.5-r-4 is sufficient. In that situation the scale where 
the top Yukawa becomes strong can be as low as ~ 100 TeV, so that the window where CT is 
active is not very big. On the other hand, the more conservative, but robust, assumption that all 
flavor violating operators are equally important at the Flavor scale requires the more constrained 

3 Of course if ds is strictly > 4 we need some other marginally relevant coupling, or a strongly relevant coupling 
which can be taken small because of an extra symmetry, to generate the weak scale by dimensional transmutation. 

4 For instance for the AS* = 2 operators contributing to KK-mixing, this amounts to assuming respectively an 
extra suppression factor ~ (V s tVdt) 2 or y s yd/y% compared to operators involving just the third family quarks. 
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pattern du — 1 & 0.2, ds > 4. That second situation corresponds to a flavor scale around 10 5 TeV, 
with the CFT describing physics in a sizeable window of scales. In view of the above, it would be 
interesting to derive, if it exists, an upper bound on ds as a function of du which was one goal 
of the present paper. Indeed in [I] a first step toward an answer was given by deriving an upper 
bound on min(ds, dx), where dx is the dimension of the triplet scalar C H x W . The method of 
[1] was however "symmetry blind" in that it could not resolving the singlet and triplet channels. 
In j2] a further refinement of the bound was obtained by working up to 18 derivatives in function 
space. The bound is a monotonically growing function of du crossing 4 at about du = 1.6. That 
result is compatible with the flavor-optimistic CT scenario. Given the clear signs of convergence 
of the symmetry blind bound [2j, and also given the remarkable success of the method in 2D 
CFTs where the bound basically tracks minimal models, we are tempted to conclude that the 
flavor-optimistic CT scenario is plausible. However the bound of [2J, if interpreted as a bound on 
ds (that is if d s < d T ), would thoroughly rule out the flavor-robust CT. 

In the present paper we have instead shown that it is possible to obtain an independent bound 
on the singlet. We have rigorously shown that the bound on ds exists and goes smoothly to 2 
as d,/, -»> 1 for (p a fundamental in SO(N) (SU(N)) and iV < 128 (N < 50). We have further 
worked it out numerically for a few small groups, and in particular for 50(4), which is the smallest 
group of phenomenological relevance. The results are not yet very strong, as seen in Table 1. In 
particular, in 50(4), ds crosses 4 already for d</, = 1.06, way within the interesting region of 
flavor-robust CT. One reason for the weakness of the bound is that our numerical method based 
on the Linear Programming algorithm does not converge fast enough when the function space is 
truncated beyond k = 4 derivatives. So our best bound just corresponds to working up to k = 4. 
One reason of the extra difficulty with respect to the symmetry blind case is that we are now 
dealing with a triple sum rule, rather than with a single one, and the complexity grows 3 2 = 9 
times faster with k. Notice that even in the symmetry blind case the bounds at low k are not very 
strong. Indeed one has the progression ~ 1.12 (k = 2) [1], d* ~ 1.18 (k = 4) (unpublished), 
d* ~ 1.35 (k = 6) p], eventually increased to d» ~ 1.6 for k = 18 [2]. Assuming a similar rate 
of improvement for the globally symmetric bound, and assuming optimistically that we could 
push the analysis to similarly high values of k (which would likely require new ideas in algorithm 
implementation), we could expect to get to ~ 1.1 -r 1.2 for the 50(4) case. This is more or less 
at the edge of interest of the flavor-robust CT. 

It is important to understand why the globally symmetric bounds are so weak, and why they 
are getting even weaker at larger N, as shown in Table 1. One could imagine two alternative 
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explanations. One, boring, possibility is that crossing symmetry for the 4-point function of just 
one operator <p a is simply not an efficient constraint in presence of global symmetry. It would 
be more interesting if, perhaps, our result is telling us something physical, namely that for larger 
global symmetry, the role of the singlet S in maintaining consistency of the theory is indeed getting 
smaller, so that it can be allowed to decouple. A very partial hint of that could be the fact, already 
emphasized in PQ, that in the O(N) model in 4 — e dimensions the anomalous dimension of the 
singlet is 0(e) while that of the symmetric traceless is 0(e/N). One way to test which of the two 
possibilities is true is to derive a twin bound on the symmetric traceless T (without imposing any 
constraints on the singlet S). If also that bound were found to be weak, the first possibility would 
be favored. If, on the other hand, the symmetric traceless bound would turn up much stronger 
than the singlet one, and perhaps comparable in strength to the general bound of [I], [2] on the 
min(dy, dg), then this would be an indication that tensor is much more important than singlet in 
maintaining the crossing symmetry. Note that, as far as we know, there's no simple reason why 
the bounds get weaker with larger N. At this stage this is just an experimental fact. 

There is yet another piece of information which, if taken into account, could change the picture 
qualitatively. In our numerical study we did not make any assumption about the stress tensor 
and symmetry current central charges. In principle our bound could become stronger under 
the condition that these central charges are bounded from above, corresponding to a perhaps 
more reasonable theory, that is one that does not contain too many degrees of freedom. Indeed, in 
connection to our main phenomenological application, the common wisdom is that constraints from 
the S parameter point towards a small EWSB sector, so that the central charges should be small. 
This is especially true of the current central charge, because the current- current spectral density 
enters directly into the spectral representation of the S parameter. Since conformal symmetry is 
broken in the IR, only the high energy tail of this density will be controlled by the CFT central 
charge. This does not allow to make this connection precise, but still large central charge seems 
to be disfavored. 

What one could do then is to study a lower bound on the central charge as a function of the 
gap in the singlet scalar sector. For the stress tensor central charge in the case without global 
symmetry, precisely such a study was performed in [1]. These studies are made possible by the 
existence of 0(1) universal bounds on the coefficients p&,i in the conformal block decomposition 
[3J, and that the coefficients p^2 and p^ t i can be related to the inverse central charges by the 
Ward identities [TT]. In [I], we found that higher gap sometimes requires a significant increase in 
the stress tensor central charge. If the current central charge is shown to have an even stronger 
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dependence on the singlet scalar gap, this may indicate a potential difficulty with the S parameter 
for CT. 

Notice that central charge studies are also interesting in their own right, without connection to 
CT. For example, can one show that an SO(N) theory with a fundamental necessarily has central 
charges larger than free theory of iV scalars? For the stress tensor central charge and the case 
without global symmetry this was shown in (in a range of near 1). 

In passing, one could imagine one day bounds like those discussed in this paper would make 
contact with the studies of IR fixed points of gauge theories performed on the lattice. Since the 
very existence of an IR fixed point implies that these theories cannot contain a singlet scalar 
with dimension below 4, our bounds could provide rigorous theoretical constraints on the lattice 
measurements of the fermion bilinear operator dimension [T2] . 

The framework laid out in this paper will likely lead to many applications beyond those men- 
tioned above. Be aware that the next crucial steps are an algorithm improvement and/or finding 
if there exists a set of questions for which the convergence is faster, so that interesting bounds can 
be obtained already at small k. This is a new field to explore! 
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